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Abstract 



The possibility of the cosmic string creation by the vacuum fluctuations of quantum fields 
in the self-consistent semiclassical theory of gravity is discussed. We use the approximate 
method for obtaining vacuum expectation value of the renormalized stress-energy tensor of 
conformally invariant quantum fields in static cylindrically symmetric spacetimes. We have 
obtained the particular solutions of Einstein equations for the different boundary conditions 
at the cylinder symmetry axis. 

There is very interesting question: can the vacuum fluctuation of quantum fields to be the 
single source of curvature or nontrivial topological structure of spacetime? One may find the 
solution of this problem in the theory of quantum gravity. But a fully satisfactory theory of 
quantum gravity does not yet exist. As a first step on the way of solution of this problem one may 
to use the self-consistent semiclassical theory of gravity. This theory sets the classical Einstein 
tensor equal to the expectation value of the stress-energy tensor operator of the quantized matter 
fields present (throughout we use units such that c = % = G = 1) 



The primary difficulty in the theory of semiclassical gravity is that the effects of the quantized 
gravitational field are ignored. A secondary difficulty of semiclassical gravity is that {T^ u ) 
depends strongly on the metric and is generally difficult to calculate. Most calculations of 
(T^) have been made on fixed classical background. These results can not be used in self- 
consistent semiclassical theory. There have been approximate calculations made for conformally 
coupled fields in any static spherically symmetric Einstein spacetime |2|. The analytical 
approximations to (T^) have been obtained by Anderson, Hiscock and Samuel || for massless 
or massive scalar fields with an arbitrary coupling £ to scalar curvature in static spherically 
symmetric spacetimes. In this work we use the Killing ansatz giving approximate vacuum 
expectation value of the stress-energy tensor (T^ u ) for conformally coupled massless field in static 
spacetimes which was obtained by Frolov and Zelnikov Q. The Killing ansatz is constructing as 
an approximate expression for the vacuum expectation value {T^ p ) using the Riemann tensor, 
the Killing vector and their covariant derivatives. The tensor (Tuv) is constructed so that the 
conservation law (TK)-„ = and the conformal trace anomaly are fulfilled. In addition a correct 
scaling transformation behaviour of (T^) is necessary. Then the following expression can be 
obtained for (T^ u ) : 






(2) 
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where g^ u is the metric tensor, £^ is the Killing vector of static spacetime, q\ and q 2 are the 
arbitrary constants of ansatz. 

The metric for the general static cylindrically symmetric spacetime is 



ds 2 



e 2F dt 2 + b 2 dz 2 + dr 2 + c z d<p 2 , 
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where F, b and c are functions of the proper distance r. 
We restrict our consideration to the special case 

F = const. 

^From this assumption it is easy to obtained the following expression for (T^ u ) 
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where the notations are used : 



e(e~ F °, 0,0,0), ^ = re CT , ^ = ^/f, 
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h (s) is the number of helisity states for the conformal massless field of a spin s. 

Four nontrivial Einstein equations for this case are given in Appendix. This equation set are 
compatible only in the case 



q 2 = 0, 2[a(F - F ) + qi] = - 
and independent equations of this set can be written in form 
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b 2 c 3bc 2 b 2 c be 2 3c 3 b 3 c be 3 I 
The solution without source at the symmetry axis must satisfy the boundary conditions: 

h" h"" 

b(r) = b + -^r 2 + -^-r 4 + 0(r 5 ) for r -» 0, 
d" 

c(r) = c' r + -^-r 3 + 0(r 5 ) for r -► 0, 
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6o,&o,Co" are arbitrary constants. These expressions can be obtained from equations (^) -(|3 
by substitution the functions b(r) and c(r) power series of r. 

The equations @ and were solved by numerical methods for the scalar field 

The results are demonstrated in figures: 
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Three cases which correspond to different boundary conditions are represented on these 
figures: 
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The asymptotic behaviour of b(r) and c(r) can be investigated analytically. From the nu- 
merical calculations, we see that for sufficiently large r, both the functions can be represented 
as a sum of two distinct components 

b(r) = b mon (r) + 5b(r); c(r) = c mon (r) + 5c(r), (14) 

where b mon (r) and c mon (r) are strictly monotone increasing, 5b(r) and 5c(r) are bounded oscil- 
lating functions. The magnitudes of these components are connected as follows: 

6h{r) - 6c{r) -c«l. (15) 



bmonij"^) C mon (r) 

Moreover, the magnitudes of the amplitudes and frequencies of 5b(r) and 5c(r) are of order 
unity: 

<$&( n >(r)~<fcW(r)~e°, (16) 

while 

bmon(r) ~ c mon (r) ~ e -1 . (17) 

In fact, one supposes that the function b mon (r) and c mcm (r) increases like r The numerical calcu- 
lations also allow us to estimate the magnitudes of the derivatives of the monotonic components: 

brnon(. r ) ~ < -rnon( r ) ~ ^ ^ - (1^) 



Having these basic estimates we can expand the equations (|9,10) keeping the leading 0(e) 
terms 

Oman oTTCt \ Z0 mon C mon J 
Cmon OTTCX \^0 mon C mon / 

The solutions of these equations are 
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where 



" 2 = iss (23) 

and Di are arbitrary constants {D\ and £>2 are complex conjugate, D%, D4, D§, Dq, D7, D$ are 
real). Note that the parameters Di are constants only in the framework of the 0(e 2 ) approx- 
imation. In fact, they may be slowly varying functions of r such that their derivatives are of 
0(e 2 ). Note also that the boundary conditions were chosen so that D3 = 0. 

In order to obtain the approximate but analytic expressions for the monotonic components 



b m on{f) and Cmon(r) more work is necessary. For this we must expand equations (32-35) up 
to order 0(e 2 ) and substitute the solutions ( ^T[ p^ ) into these equations. Take into account 
that these equations contain terms of three different types: the terms which consist only of the 
monotonic functions b mon {r) and c mon (r); the terms which contain one oscillating function 5b(r) 
or 5c{l) (or its derivatives); and the terms which contain a product of two oscillating functions 
8b(r) or 5c(l) (or its derivatives). After substituting the solutions (^J2^) into (|3^-^) one can 
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check that the these equations consist of two independent parts: the monotonic one and the 
oscillating one. Both these parts must be equated to zero separately. Equating the monotonic 
part to zero we obtain the equations for b mon (r) and c mon (r): 
the 00 component 
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the zz component 
the rr component 
the ipcp component 
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And the oscillating part of these equations gives the correction of the order e 2 for 5b(r) and 
5c(r). The solutions of equations ( 24 -p7|) are 
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+ fc 4 + 0(e), 
! +0(6 3 ) : 



(28) 
(29) 
(30) 



where fcj are constants. 

Finally, the asymptotic form of metric is 
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The asymptotic form of the angular part of these solutions coincides with the angular part 
of cosmic string. However, in the case of string the conic part of spacetime is 2-dimensional 
Minkowski spacetime minus a wedge and in the considered cases the conic part of spacetime is 
2-dimensional Minkowski spacetime plus a wedge. The value of the cone excess angle is of order 
unit and it is determined by the boundary conditions at the symmetry axis. 
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Appendix 

The equation G§ = 8vr(T ) give 

b l + d ! + ^ - 8vra _ ^ + _ 1^ _ 4c "' 6 ' 

6 c 6c \ 3i 3c 36 2 36c 36c 

2c"y 46" 2 6"c" 4c" 2 26"6' 2 46"6V 26V 2 2c"6' 2 
+ lk 2_ + ~W ~ 2 ~6c" + 3c 2 " " IP" + IPF" + "36?" + l^" 

4c // b / c / 2cV 2 26 /3 c' 26 /2 c /2 26V 3 \ 
+ 36c 2 3c 3 36 3 c + 36 2 c 2 36c 3 ) ' 

The equation G\ = 8vr(T 2 2 ) give 

d' _ ( 46"" 2c"" 8b"'b' 8b"'c' 2c"'b' 2d" d b" 2 26" c" 4c" 

_ _ s^a ^--3^- + -^r + §pr + 3^2- + 2 "5T " ~ ^ 

86"6' 2 26"6 , c' 46" c' 2 4c"6' 2 8c"6'c' 2d'd 2 46 ,3 c' 26 ,2 c' 2 46V 3 \ 
W + ~WT + "36c 2 3W + "36c 2 " + ~lk^~ + ~Wc~ ~ 'Wd 2 ~ ) ' 

The equation G r r = 8vr(T r r ) give 

b'd _ ( 4b"'b' 2b'" d 2d"b' 4d"d 2b" 2 2b" d' 2d' 2 4b"b' 2 
— _ 8ira y-^yr + + "3^ 3^2" + ^ " + + 

b"b'd 2b" d 2 2d'b' 2 d'b'd Ad'd 2 2b' 3 c 2b' 2 d 2 2b'd 3 \ 
~ 2 l?c~ + "36c 2 " + ~Wc~ ~ 2 lbc T + Ik 3 " ~ 3637 + ~WJ ~ "36c 3 " ) ' 

The equation G£ = &ir(Tg) give 

b" _ (2b"" 4c"" 2b'" b' 2b'" d 8d"b' 8d"d 46" 2 26"c" c" 
~b ~ \Tb W + ~3bc 36c~~ + Ik 2 "" !ibc~ + 2 ? 

26"6 /2 8b"b'd 46"c /2 4c" 6 /2 2c"6 / c / 8c"c /2 46^V 26 /2 c /2 46 / c /3 \ 
+ ^63~ + ~Wc 36c 2 " + ~Wc~ + Ike 2 Jk 3 "" ~ ~Wc~ ~ "36V" + "36c 3 " ) ' 
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